Hawking radiation in dispersive theories, the two regimes 
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We compute the black hole radiation spectrum in the presence of high frequency dispersion in a 
large set of situations. In all cases, the spectrum diverges like the inverse of the Killing frequency. 
When studying the low frequency spectrum we find only two regimes: an adiabatic one where the 
corrections with respect to the standard temperature are small, and an abrupt one regulated by 
dispersion which is well described using step functions. The transition from one regime to the 
other is governed by a single parameter which also governs the net redshift undergone by dispersive 
modes. These results can be used to characterize the quasi-particles spectrum of recent and future 
experiments aiming to detect the analogue Hawking radiation. They also apply to theories of 
quantum gravity that violate Lorentz invariance. 
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I. INTRODUCTION 

For long wavelengths, the propagation of sound waves 
in a moving fluid is analogous to that of light in a curved 
metric [1]. However, Hawking radiation issues from very 
short wavelength modes [2-5]. Therefore dispersion ef- 
fects must be taken into account when computing the 
spectrum emitted by an acoustic black hole. When the 
dispersive scale £ (the healing length in a Bosc conden- 
sate) is much smaller than the surface gravity scale 1/k 
which fixes the Hawking temperature Th = k/2tt (in 
units where h = &b = c = 1), it is now well under- 
stood that the spectrum is robust [6-17]. However the 
precise extension of this regime, and the spectral proper- 
ties outside it are still poorly understood. 

The aim of this paper is to remedy this situation. Do- 
ing so we shall complete Ref. [15] where we studied flow 
profiles with a local perturbation that strongly varies 
close to the horizon. In that paper we found two regimes: 
as long as the perturbation scale is larger than a certain 
critical length, the temperature is the standard one fixed 
by k. For variations on shorter distances instead, the 
temperature is given by a spatial average over that criti- 
cal length. Interestingly the latter is not simply given by 
the dispersive length £. Rather it scales as £ 2 / 3 k -1 / 3 for 
quartic dispersion. 

To complete that analysis, we study the spectrum in 
simpler flows characterized by less parameters. We ob- 
serve two well denned regimes. Firstly, when the flow is 
smooth and dispersion negligible, the spectrum slightly 
differs from the relativistic one. Secondly, there is an 



opposite regime of strong dispersion and steep horizon 
which is well described by using step- like profiles [18, 19]. 
With numerical techniques [14-16], we investigate the in- 
termediate situations and show that there are no other 
regimes than these two, in agreement with [20]. We also 
show that the transition from one regime to other is gov- 
erned by a single parameter which scales in the same way 
as the above mentioned critical distance. Interestingly, 
this parameter also fixes the (finite) redshift undergone 
by dispersive modes when scattered on a sonic horizon. 
For definiteness, we consider the phonon spectra in the 
context of atomic Bose condensates. Yet, our analysis ap- 
plies to other superluminal dispersion relations, as well 
as subluminal ones, in virtue of the correspondence [17] 
between the various types of dispersion relations. It can 
therefore be used to explore the consequences on black 
hole thermodynamics [21] that arise from violations of 
Lorentz invariance found in certain theories of quantum 
gravity [22-25]. 
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The paper is organized as follows. We first present the 
relevant parameters that govern the mode scattering on a 
sonic horizon. We then compute the low frequency limit 
of the emitted fluxes in step-like flow profiles. Finally, by 
numerically scanning a wide parameter space, we char- 
acterize the validity domains of both the Hawking and 
the step-like regimes. In Appendixes we present the ba- 
sic concepts governing phonon fluxes, as well as the tools 
needed to compute these fluxes in step-like flows. We 
also show that our main results remain unchanged when 
considering a larger class of flows. 
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II. RELEVANT CONCEPTS 

A. Number of e- folds and critical frequency w ma x 

We consider elongated condensates stationary flow- 
ing along the longitudinal direction x, and we assume 
that the transverse dimensions are small enough that 
phonon excitations are longitudinal. In the hydrody- 
namic regime, phonon propagation is governed by the 
d'Alcmbert equation in the metric [1, 26] 



ds 2 = -c 2 dt 2 + {dx~v dtf 



(1) 



where the flow velocity v and the speed of sound c only 
depend on x. We assume that the condensate flows from 
right to left (v < 0), and that the sonic horizon (where 
w = c + v crosses 0) is located at x = 0. When the 
gradient 



a, 



W\ X= Q 



(2) 



is positive, one obtains an analogue black horizon. In- 
deed, in the near horizon region (NHR) where w ~ 
kx is a good approximation, the null geodesies are 
x = xoe K ^~ to \ Correspondingly, the wave vectors 
obey [5, 27] 



k(t) = fc e- K( *-* >. 



(3) 



These equations are the signatures of a (Killing) horizon. 
Notice that this exponential decay law is also found in an 
inflationary dc Sitter space where the Hubble parameter 
H plays the role of k. In that case, the number of e-folds 
TV = H(t — to) is used to quantify the duration of the 
inflationary period. The origin of this instructive corre- 
spondence is simple: in the NHR, the metric of Eq. (1) 
coincides with that of the flat sections of de Sitter in 
Lemaitre coordinates, for details see Sec. 3.3 in [28]. 

When ignoring short distance dispersion, Eq. (3) im- 
plies that n^, the spectrum of upstream phonons spon- 
taneously emitted from the horizon is Planckian with 
Hawking temperature Th = k/2tt, in units where Ti = 
ks = 1- When introducing dispersion, the spectrum n u 
is more complicated. To characterize its properties, we 
define the temperature function T M by 



1 



cxp(w/T w ) - 1 ' 



(4) 



As such, T u is simply another way to express n u . It is 
particularly useful here because n u always diverges like 
1/ui for oj —> [14, 15, 18]. Hence the low frequency part 
of the spectrum is governed by the effective temperature 



T eS = lim T u 



lim to n u 

uj->0 



(5) 



We define the Hawking regime by situations where T e ff 
differs from Th by less than 10%. We aim to determine 
the extension of this regime when introducing short dis- 
tance dispersion. 



In a stationary flowing atomic Bose condensate, the 
dispersion relation of phonons is given by 



(uj-vk) 2 = n 2 {k) 
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(6) 



where uj is the conserved frequency idt , k the wave vector 
—id x , m the atom mass, and A the frequency associated 
with the healing length £ = Ti/ y/2mc by A = \/2c/£. This 
dispersion relation directly follows from the Bogoliubov- 
de Gennes equation [29], see Eq. (A6) in Appendix A, 
where more details can be found. Because of dispersion, 
phonons wave packets no longer follow null geodesies in 
the NHR [7, 12]. However, as long as the packet prop- 
agates in that region, Eq. (3) applies, irrespectively of 
the value of A/k, and in fact irrespectively of the dis- 
persion relation f2(fc). What is essential here is that the 
trajectories stay in the NHR only for a finite time. As a 
result, unlike what is found in de Sitter space, the accu- 
mulated red shift effect saturates at a finite value (when 
the asymptotic values of c and v arc finite, which is the 
case we shall consider). To quantify this saturation, we 
use 



„(2),as' 



(7) 



where ku 33 and fc^ 2 ^' as are the asymptotic values of the 
in-going and out-going wave numbers, solutions of Eq. (6) 
with f2 < 0, see Appendix A. These roots describe the 
negative frequency partner which is trapped in the super- 
sonic region, on the left of the horizon. As we shall see, 
the number of e-folds N Ml or its exponential r u , offers a 
simple description of the low frequency spectrum outside 
the Hawking regime. It is thus important to determine 
how it depends on the dispersive scale A and on the flow 
profile parameters. 

To reduce the number of these parameters, we impose 



c(x) = c + (1 - q)w{x), 
v(x) = -c + qw(x), 



(8) 



where Co is the sound speed at the horizon. These equa- 
tions specify how w = c + v is shared between c and v. 
In the forthcoming numerical simulations, we shall use 
q = 0.3, which is a representative value of the generic 
case [14], and w given by 1 



w(x) _ / KX 

— - — - = D tanh 

c \c D 



(9) 



The important quantity here introduced is D. It fixes the 
asymptotic value of w, but, more importantly, it also fixes 



1 In Appendix C, we consider more general profiles and observe 
that the main results are not significantly modified. 



3 



the extension of the NHR, i.e., the domain where w ~ 
kx. This latter fact renders D a very relevant spectral 
parameter in the presence of dispersion [15, 17]. 

To compute r u of Eq. (7), we need to analyze the two 
roots of Eq. (6) with Q < 0. In the supersonic region, 
these roots merge in a point that corresponds to x tp (ui), 
the turning point of the trajectory at fixed u> [17]. Its 
location is fixed by solving v gr = 0, where the group 
velocity (in the horizon rest frame) is 



v S r = dkLo = 9^(0 + vk). 
Using Eq. (6), for uj <C A, one finds [8] 



CO 



OA 2/3 

A/ 



(10) 



(11) 



For cj /A — > 0, x tp coincides with the horizon at x = 0, 
as expected from the behavior of the trajectories in the 
absence of dispersion. When u increases, the turning 
point recedes to the left. We define the critical frequency 
w max by the value of u> where x tp (u;) reaches — oo. In the 
velocity flows of Eq. (9), one has w max = A f(D, q) where 
/ is a rather complicated function. For small values of 
D, which will be a relevant regime in the sequel, one gets 



A 







D 






or 


H 




H) 


+ 0(D 2 ) 






' ~2 







(12) 



Using the asymptotic values of v and c, we can compute 
r u of Eq. (7), and relate it to w max . After a lengthy 
computation we obtain 



B. Fluxes for step-like profiles 

When ignoring grcybody factors, the spectral proper- 
ties of riuj are analytically known in the Hawking regime 
which is found in the dispersionless limit (A — > oo, 
fixed ft). In the opposite limit of infinite surface grav- 
ity (ft — > oo, fixed A), the spectral properties can be 
analytically computed [18-20]. In particular, the low fre- 
quency effective temperature T e g of Eq. (5) is well de- 
fined. We call it T s t C p and we refer to this case as the 
step-like regime. 

Before determining the extension of this domain, we 
generalize the treatment of the above references so as to 
work with condensate flows where both v and c have a 
jump at the sonic horizon. We focus on the low-frequency 
behavior of n w , the spectrum emitted in the subsonic re- 
gion, and n^j, the spectrum of left movers with respect 
to the atoms, see Eq. (A10). In the step-like regime, un- 
like what is found in the Hawking regime, the spectra of 
both right and left movers are predicted by the same tech- 
niques and should be treated on equal footing. In fact, 
whereas the spectrum of left-going particles is not well 
defined for gravitational black holes as there is no asymp- 
totic region inside the horizon, it is perfectly well defined 
for acoustic black holes. Moreover the mode mixing of 
left and right movers has been observed in gravity waves 
experiments [30] and could be observed in the future in 
Bose-Einstein condensates. In Appendixes A and B we 
give a summary of the concepts needed to compute these 
two occupation numbers. Here we present only the main 
results. 

In the low frequency limit, we obtain 



3^ 



D 
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3V3 
0(D 2 ) 
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1 + T 
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U) 



(13) 



For frequencies co/uj max <C 1 and small values of D, quite 
remarkably, r w depends only on u) max /u). 

At this point we remind the reader that in dispersive 
theories it was observed that w max governs the first order 

deviations from the Hawking spectrum [13 16]. Since 

(i) 

the two negative roots ku , associated with the negative 
norm mode, no longer exist for w > o> max , w max acts as a 
cutoff above which there is no radiation. However, it is 
more subtle to understand why w max also determines the 
low frequency (cj -C w max ) properties of the spectrum. A 
detailed analysis [17] of the connection formula governing 
the scattering on the sonic horizon confirms that the first 
deviations from the Hawking spectrum are governed by 
the particular combination of A and D entering in w max . 
From Eq. (13) we learn that the same combination fixes 
the total redshift when D is small. As a result, we shall 
see that the transition from the Hawking regime to the 
other one can be meaningfully characterized by r u for 
lu = T cS of Eq. (5). 
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Using Eqs. (5) and (12), the effective temperature reads 
3V3 



step 



i-U- g )i>- 



O D 2 



(15) 



This simple equation establishes the relevance of w roax 
in the step-like regime. In particular it establishes that 
n^j cx D 3 / 2 as D — > 0. This non-trivial result is relevant 
for experiments in optical fibers and glasses where D < 
10 -3 . It is also interesting to express n^ tGp in terms of r w 
of Eq. (7): 



step 



^[l- qD + 0(D 2 )] 



(16) 



Hence, for small values of D, using Eq. (B18), we get 
another simple equation: 



step 



4c 



(17) 



where v u is the fluid velocity in the upstream region. We 
notice that n^ tep factorizes into r w and v u , depending on 
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quantities which are respectively defined in the super- 
sonic and subsonic region only. 

Operating in a similar way, we obtain the number of 
left-going particles 



7>, Step 



q 

D 
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A 



D 7/2 



q+ 2 



0(D 2 ) 



+ Cst. (18) 



In analogy with Eq. (5), we define T e " ff the effective tem- 
perature of left movers. In the step-like regime, we get 



rpV 

step 



3^3 



D 2 



q -2 



0(D 4 



(19) 



Note that T s " cp vanishes quadratically when q = 1/2. 
This is reminiscent of what was numerically noticed 
in [14] for continuous profiles. In that case it was found 
that the mixing between right- and left-going modes was 
minimum for this value. For a step-like profile, the first 
non- vanishing term is 



v, step 



D 11 ' 2 
1152^2 



1 



D 
~2 



A 



O [ D 13 ' 2 



(20) 



Finally, since both and n w diverge in the same way 
for uj — > 0, it is convenient to compute their ratio 



Rcff 



lim ^ 



(21) 



For step- like profiles, using Eqs. (14) and (20), one gets 



stop 



4 [q 



2 1 ° 2 



O (D 



(22) 



From this ratio, we learn that is highly suppressed 
in optical fibers since D <C 1, and that it is not directly 
related to r u of Eq. (7), unlike what is found for the right 
movers in Eq. (16). In that case indeed, for D <C 1, both 
T s tep and r u have the same dependence on A, D and q. 

In the next section, we determine the validity domain 
of these formulas. It is rather clear that the above ex- 
pressions should be valid in the limit of infinite surface 
gravity: k/A — > 00 at fixed D. Rather surprisingly, we 
shall see that they are also valid for very small k/A when 
D 0. 



III. THE EXTENSION OF THE TWO REGIMES 

A. The Hawking regime 

We start with the extension of the Hawking regime. 
This has been already studied and it is rather well under- 
stood. Using numerical techniques [14], it was found that 
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log(A/ K) 

FIG. 1: Teff/TH as a function of A/k and uw x /k; in logarith- 
mic scales. The parameter q of Eq. (8) is 0.3. Solid lines are 
loeii of constant T c g/Tn- One clearly sees that the departures 
from the Hawking regime are essentially governed by u max /K. 
Roughly speaking, one leaves the Hawking regime (i.e., T e g 
and Th differ by 10%) for oj max /K < 0.4, irrespectively of the 
value of A/k. N.B. The upper white area is due to the fact 
that we restrict our analysis to D < 1, whereas the lower 
one is due to our incapacity of numerically obtaining n u for 
D < 0.02. 



the first order deviations from the Hawking spectrum are 
mainly controlled by w max /K. They thus scale as AD 3 / 2 , 
see Eq. (12). In fact, it was found that these deviations 
hardly depend on A/k when working at fixed w max //«. 
Using the same code, we now clarify these observations 
and also investigate the spectral properties further away 
from the Hawking regime. To these ends, we plot in Fig. 1 
the ratio T e g /Th in the two-dimensional parameter space 
defined by w max //« and A/k. We use logarithmic scales 
to see more easily the scaling properties. The value of 
the parameter q has been chosen to be 0.3, to avoid the 
peculiar cases ofg=l/2or(jr=l/6, where some contri- 
bution cancels, as can be seen from the expressions given 
in Sec. II B. 

For w max /«; > 1, we see that T e g/Tn — > 1 in a manner 
that hardly depends on A/k. The region o; max / k > 1 
therefore characterizes the situations where Hawking's 
result is robust. For w max / ' k <C 1, the effective tempera- 
ture T e ff, strongly differs from Th- However it still hardly 
depends on the dispersive scale A/k. Therefore U) max /K is 
the most relevant spectral parameter in the entire plane. 
This important lesson is corroborated by the results pre- 
sented in the following subsection. 
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FIG. 2: The temperature ratios T e g/T s tc P (left panel) and T^ s /T^ tcp (right panel) as a function of A/k and t<j max /K for the 
same situations as in Fig. 1, where T s t ep and T s " ep are obtained from Eqs. (B16) and (B17). The departures from the step-like 
regime are essentially governed by Lo ma _ x /n. Roughly speaking, one leaves this regime for c^max/^ 0.1, irrespectively of the 
value of A/k. 



B. The step-like regime 

In a manner strictly similar to what we have just done 
for the Hawking regime, we study the validity domain of 
the expressions obtained in Sec. II B. In Fig. 2 we plot 
T cS /T stcp (left panel) and T^/T£ ep (right panel) for the 
same situations as in Fig. 1. The first lesson from Fig. 2 
is that these two ratios behave in a very similar manner. 
Secondly, when u; max / k < 0.1, we see that the effective 
temperatures T ff and T£n obtained with the regular pro- 
files of Eq. (9) are in excellent agreement with T step and 
T s ^ ep obtained in Sec. II B. In addition, we learn that this 
agreement is reached in a manner that hardly depends on 
A/n. Therefore, quite remarkably, the step-like regime is 
found even when scale separation is achieved, i.e., for 
A/k ^> 1, when D is small enough that w max /K < 0.1. 

To complete our analysis, in Fig. 3 we plot the ratio 
R e g of Eq. (21). We note that the lines of constant R e g 
almost coincide with the lines of constant D (slope 2/3, 
dashed lines). Hence R c g scales as D 2 . This scaling 
is in perfect agreement with the step-like prediction of 
Eq. (22). This is a further confirmation of the applicabil- 
ity of step-like techniques in a wide region of parameters. 
However, unlike T e ff and T^ ff that are well approximated 
by using a step-like profile only for w max /«; < 0.1, their 
ratio behaves as D 2 in the whole plane. This means that 
i? e ff scales in the same way in the Hawking and in the 
step-like regime. 

We also emphasize that R c r is maximal for large values 
of D. This is also rather unexpected as it was believed 



that the coupling between right and left sectors should be 
a sub-dominant effect for smooth profiles. When adding 
high frequency dispersion to a 2D conformally invariant 
masslcss field [6, 7], this is indeed the case [13]. However 
when using the quasi-particle wave equation of some con- 
densed matter model, there is no reason to believe that 
this should still be the case. (From Figs. 11 and 12 of 
[14] it is clear that in Bose-Einstein condensates, the 
spectrum of left movers is generically large.) 



IV. SUMMARY AND INTERPRETATION 

To summarize the results, we schematically illustrate 
in Fig. 4 the validity regions of the two regimes. The 
region where the temperature T c ff differs from Th by less 
than 10% is shaded in dark gray, whereas that approxi- 
matively governed by a step-like profile is shaded in light 
gray. We re-emphasize that the expressions obtained us- 
ing a step-like profile provide very good approximations 
not only where k/A is large (right lower corner), but in 
general when w max / k is small, independently of the value 
of k/A. 

In Fig. 4, we see that the white area from 0.1 < 
u max /K < 0.4 is not covered by either approximation. To 
characterize the effective temperature in this transitory 
regime, we use the average temperature T av of [15]. It is 
obtained by taking the spatial average of the gradient of 
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log(A/ K) 

FIG. 3: Tieff defined in Eq. (21) as a function of A/k and 
Wmai/c for the same situations as in Fig. 1. Solid lines give 
the various values of this ratio, whereas the dashed straight 
lines with slope 2/3 are locii of constant D, see Eq. (12). 
Clearly, i? c ff depends only on D. 
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FIG. 5: Tefj/Tav as a function of A/k and o; max /fi: in loga- 
rithmic scale for the same profiles as in Fig. 1, where T av is 
given by Eqs. (23) and (24). Solid lines are locii of constant 
T e ff/TB. V . The other lines are as in Fig. 4. 
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FIG. 4: In dark gray we denote the validity region (defined 
by 10% difference in temperature) of the Hawking regime, and 
in light gray that of the step- like regime. The thick horizontal 
line is the threshold value of Eq. (26). The dashed straight 
lines are locii of constant D. It is clear that cj max /K governs 
the extension of the two domains, and the transition between 
them. 



w = c + v over a given length 

[ d t /2 dw _ 2w(d € /2) 



dx — — 

dx 



-d € /2 « x «£ 
The critical length was found to be 



d£ = dco (kA 2 



-1/3 



(23) 



(24) 



where d is approximatively equal to 2 2 / 3 . When the ex- 
tension of the NHR where w ~ nx is larger than d^ 1 T av 
agrees with the Hawking temperature Tjj. 

Transposing these notions to the present context, the 
boundary of the region where the Hawking regime ceases 
to be valid should thus be given by the condition that 
the extension of the NHR of Eq. (9) equals d^: 



2^ = d,. 

K 



(25) 



Using Eqs. (24) and (12), this condition fixes the thresh- 
old value of Wmax/ft to be 



threshold 



3V3' 



(26) 



In Fig. 4, the black solid line represents this value. In 
accord with the analysis of [15, 17], it precisely lies where 
the Hawking regime ceases to be valid. 

To be more precise and more quantitative, we now 
compare T av with the effective temperature T e g. Their 
ratio is plotted in Fig. 5. As expected, T av is very close to 
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FIG. 6: The number of e- folds of Eq. (7) evaluated for uj — T cfi 
as a function of A/k and w ma x/K for the same profiles as in 
Fig. 1. The dotted lines bound the regions of validity (with 
10% errors) of the Hawking regime (upper line) and the step- 
like regime (lower line). The thick line is the threshold value 
of Eq. (26). 



Secondly, we point out that our analysis can be easily 
used to characterize the outcome of recent and future ex- 
perimental investigations aiming to detect the analogue 
Hawking effect. Indeed the knowledge of A/k and w max /fv 
suffices to determine in which regime one will find the 
(low frequency part of the) quasi-particle spectrum. For 
instance, a step-like analysis should be appropriate for 
experiments in optical fibers or glasses since these are 
characterized by a small parameter D, typically of order 
< 1(T 3 [31, 32], in spite of the fact that A/k > 1. Thus, 
in these cases, a precise knowledge of the surface gravity 
is not required to predict the spectrum. On the contrary, 
for experiments in Bose-Einstein condensates [33] or in 
water [30], the spectrum will be essentially governed by 
the surface gravity because D is of the order of unity, and 
A/k > 1. 

Thirdly, our analysis can be used to characterize the 
Hawking radiation in theories of quantum gravity [22-25] 
which contain violations of Lorentz invariance at very 
high energies. In particular, Fig. 1 indicates that the 
black hole temperature should no longer be given by the 
surface gravity when n/A becomes larger than 1. Un- 
derstanding the possible consequences of this on black 
hole thermodynamics [21] is an interesting and challeng- 
ing question. 
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T e g when D and A/k are large. When D <C 1 , T av is a bit 
less accurate in reproducing T e g than Th, even if u^ m ax/ K 
is large. However, T av gives a rather good approximation 
of T e ff in a much broader region than Tfj. Within an 
error of 10%, the region covered by T av extends down to 
log 10 u> max / K Rj — 1. The average temperature T av there- 
fore correctly describes T c g in the transition region where 
both Th and T st0 p fail to do so. 

Now, what is the physical meaning of Eq. (26)? Using 
Eq. (13), we see that w max is proportional to wr u . To be 
able to define a threshold value for the number of e-folds 
N u = lnr w , we need to pick a frequency u. In the sequel, 
we work with oj = T B g since T e g characterizes the actual 
spectrum n u for uj -> 0. When w max < w max | threshold, we 
see in Fig. 6 that N u —T e{i is almost constant and slightly 
larger than 1. It has to be noticed that it does not go 
to as D — > when working with uj = T c ff . When cj max 
is above the threshold, Nt o{{ increases exponentially as a 
function w max /K. We therefore conclude that the number 
of e-folds evaluated at the typical frequency T c g provides 
an extremely simple way to characterize the validity do- 
mains of the two regimes. That is, the Hawking regime 
starts above Nt cH ~ 2, whereas the step- like analysis 
provides reliable results when Nt b!! is smaller than 1.4. 

To conclude this paper we first underline that our re- 
sults, which have been obtained using the superluminal 
quartic dispersion of Eq. (6), also apply to other types of 
superluminal relations, and to subluminal ones, in virtue 
of the correspondence between these cases [17]. 
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Appendix A: Mode analysis and scattering matrix 

We present the main concepts that characterize the 
scattering of phonons in atomic Bosc condensates in sta- 
tionary flows containing one sonic horizon. We follow the 
treatment of Ref. [14] where more details can be found. 

In non homogeneous condensates, phonon elementary 
excitations are appropriately described by relative den- 
sity perturbations: 

<F = *o(l+0), (Al) 

where ^ is the atom field operator, and VPrj the con- 
densate. The dynamics of <j> is determined by the 
Bogoliubov-de Gennes equation [29]: 

id t (j> = [T p - ivd x + mc 2 ] $ + mc 2 $\ (A2) 

where Ti = 1, c is the x-dependent speed of sound 
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Po(x) = \^>q(x)\ 2 is the mean density of condensed atoms, 
and g(x) is the effective coupling constant among atoms. 
These functions play no role in the sequel: only c(x) and 
v(x) matter. The kinetic operator that acts on <j) is 



-—vd x -d x . 
Zm v 



(A4) 



In stationary situations, 
quency eigenmodes 



can be expanded in fre- 



t <Mx)+e+ K ^(x)t 



= J duj 

= /^E [^'CM^ + e + ^(vZ(x))*aZ 1 ] , (A5) 

J a 

where the discrete sum over a takes into account the 
number of modes at fixed uj. Inserting Eq. (A5) in 
Eq. (A2) yields the following system 



[(w + ivd x ) -T p - mc 2 ] 4>Z( X ) = m c 2 ip"(x), 
[-(w + ivd x ) -T p - mc 2 ] (p"(x) = mc 2 (f)"(x). 



(A6) 



In backgrounds containing one sonic horizon, there are 
two or three modes depending on whether uj is larger 
or smaller than w max of Eq. (12). More precisely the uj 
component of the field operator reads 



\ * 2.1 



(A7) 



For uj > 0, the modes 0™ and <f>^ have positive norm, 
they describe, respectively, propagating right-it and left- 
v moving waves with respect to the condensed atoms and 
they are associated with the usual real roots of Eq. (6), 
see Fig. 7. On the contrary, (<p_ u )* has a negative norm. 
It is associated with the extra roots that exist in super- 
sonic flows for uj < uj max , and describes phonons that 
are trapped inside the horizon (x < 0) in the supersonic 
region. 

In our infinite condensates, the above modes become 
superpositions of plane waves both in the left and right 
asymptotic regions of Eq. ( A7) . We can thus construct 
without ambiguity the in- and out-mode bases. The in 
(out) modes are such that each of them contains only one 
asymptotic branch with group velocity directed towards 
(away from) the horizon. For uj > w max , there exist only 
two positive norm modes, and therefore in and out modes 
are linearly related by a trivial (elastic) 2x2 transforma- 
tion. Instead, when uj < w max , the in and out modes mix 
with each other in a nontrivial way by a 3 x 3 Bogoliubov 
transformation [14] 

C' in = «-C' out + + ^C out , 

(^-J* = ^C'° ut + a- u (tp™*)' + ^C° ut , (A8) 

The standard mode normalization yields relations such 
as (from the first equation) 



\a u 



(A9) 



When the initial state is vacuum, the final mean occupa- 
tion numbers are given by 



1/3. 



2 

uj | 1 



n-u, = n. 



(A10) 



The third equation, which gives the number of nega- 
tive frequency phonons, follows from energy conservation 
which imposes that each pair spontaneously produced 
should have no energy. In the standard analysis without 
dispersion, is Planckian and describes the Hawking 
effect [5]. In that case, for massless conformally invari- 
ant fields, the number ri^ of left moving phonons iden- 
tically vanishes because the coupling between u and v 
modes is zero. Adding high frequency dispersion to that 
case, one finds -c n u . Hence the scattering of this 
dispersive field on a sonic horizon basically consists of a 
more general and slightly modified version of the stan- 
dard Hawking effect, as long as w max > k. For larger 
values of k, and for non-conformally invariant fields, de- 
viations with respect to the rclativistic case can be large, 
but the structure of Eq. (A8) and the meaning of its co- 
efficients remain unchanged [14, 18, 19]. 



Appendix B: Step-like profiles 

We extend the computation of the Bogoliubov coeffi- 
cients given in [18, 19] to step-like profiles where the ve- 
locity is different on the upstream and downstream sides 
of the horizon. That is we consider profiles given by 



v(x) 
c{x) 



v u 9(x) +v d 0(-x), 
c u 6(x) + c d 0(-x). 



(Bl) 
(B2) 



For these profiles, Eq. (A6) can be solved separately in 
the upstream (x > 0) and downstream (x < 0) regions 
by Fourier transform: 




-iut+ikiX 



duj ^/Wpl 

dk g—iwt+ikix 

duj y/Wpi 



(B3) 
(B4) 



where the index i stands for u (upstream) when x > and 
d (downstream) when x < 0. (For the sake of conciseness 
we have omitted the indexes a and uj). The normalization 
factor has been chosen such that 

u 2 - = 1, (B5) 

where Ui, Vi are chosen real. On each side Eq. (A6) gives 



k 2 

(uj - Vi h) - — l - mc 2 

Ira 

k 2 

-(uj - v,i hi) - - 1 mc 2 

2m 



mc 2 Vi 



mc 2 Ui. 



(B6) 
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0,/k Vl/n 




FIG. 7: Graphical solution of the dispersion relation (6) for subsonic (left panel) and supersonic flow (right panel). Solid 
curves: ±fi(fc)/«. In subsonic flows, for to > the dashed line (u> — vk)/n crosses twice the dispersion relation with Q > 0. In 
supersonic flows, two extra roots exist in the left lower quadrant for u < o; max (dashed line). For u — cj ma x the solid line is 
tangent to fi(fc), and for to > (xJ max (dotted line) these real roots no longer exist. 



Since the dispersion relation is quartic, we obtain four 
roots k(uj). Hence, for a given frequency ui, the general 
solution can be written on either side as superposition of 
four momentum cigenmodcs. Notice that k are all real in 
the downstream supersonic region, whereas two are real 
and the other two complex conjugated in the upstream 
region. 

In order to obtain globally defined solutions we need 
to match the two sets of modes. Firstly, the solutions 
4>u,i{x) and (fiui,i(x) must be continuous at the horizon, 
that is 



(B7) 



where the superscript ± denotes the limit of the func- 
tion for x — > 0^, respectively. Secondly, the matching 
conditions for the derivatives of <fi and ip can be obtained 
by integrating the system (A6) in a small neighborhood 
[—£, e] of the horizon and taking the limit for e — > 0. To 
characterize this limit, it is useful to write T p of Eq. (A4) 
as 

1 dl + J-[d x \og(v(x))]d x . (B8) 



2m 2m 
For our profiles, the second term gives 

d x \og{v{x)) = d x [logK) 9{x) + logM 6(-x)} 

= log(v u ) 6(x) - \og(v d ) 6(-x) = log (^j S(x). (B9) 

The non-trivial terms of the first equation of Eq. (A6) 
are 

i f dxvd x 6-\ / dxdtd) 

J-e 2mi_ e xV 

--Llog(^i) / E dx5{x)d x <j> = 0. (BIO) 
2m \v d J J_ e 

Since </> is continuous at x = 0, its derivative can have 
at most a finite jump at x = 0. Therefore, the integrand 
in the first integral has a finite jump at x = 0, such 



that, when taking the limit e — > 0, this term vanishes. 
The second integral can be directly computed because 
its integrand is the derivative of d x 4>. It gives the jump 
in the first derivative of the field d>: 



(Bll) 



dxdl4> = d x <\>X - d x cj> d 
Finally, the last integral gives 



log^^-^ J dx5(x)d x (j) 



Taking into account these last equations, the matching 
condition of the first derivative is 



l--log(- 

2 & \v d 



1 + - log — 

2 \v d 



Or 



(B13) 

From the second equation of ( A6) , one finds an identical 
condition for d x tp. In brief, we get 



where 



V 



_ 2 - log(v d /v u ) 
2 - \og(v u /v d ) ' 



(B14) 



(B15) 



Using Eqs. (B7) and (B14), for ui < w max , one can deter- 
mine the three globally defined solutions that are asymp- 
totically bounded and that appear in Eq. (A7). From 
them we can extract the nine coefficients of Eq. (A8). 
The computation is rather tedious, but completely anal- 
ogous to what was performed in [18, 19]. Below we only 
give the dominant term in the limit lu — > of two Bo- 
goliubov coefficients f3 u and B u which determine n w and 
n v ,. We obtain 
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A 

x 

CO 

A 



lv_d jc 

v 



U I u u 



+ c u (v u - v d /ri) 



VdV u - V c d 



U 



uj 2c 0y /cd~ V Cu-v-u [cl - c 2 d + c u {v u - Vd/rj)] 



Cu/V + v d 



Cd — — c5 1 



CdV + Vu 

Cd + Vd 



2 _ „2 



Cd 



(B16) 



(B17) 



Using our parameterization of the velocity profile given 
in Eqs. (8) and (9) 



c u = c [1 + (1 - q)D] , v u 
Cd = c [1 - (1 - q)D] , v d 



-co(l-gD), (B18) 
-CD(l + gD), (B19) 



in the limit of small D, we obtain Eqs. (14) and (20). 



Appendix C: Asymmetric profiles 

In this Appendix we show that the results obtained in 
the main body of the text are still valid when considering 
flow profiles that generalize Eq. (9): 



w a (x) 



D 



Co 



/ KX 

a + tanh 

\c D 



(Cl) 



The parameter a fixes the asymmetry between the sub 
and supersonic flows with respect to the sonic horizon. 
This horizon exists only if \a\ < 1, and it is located at 



ih = arctanh(a). 

K 



(C2) 



The Hawking temperature corresponding to the velocity 
profile of Eq. (Cl) is Th = K a /2ir, where 



K a — d x W a \x— x 



K [1 



(C3) 



In [16], we showed that when \a\ < 1, the spectrum 
remains Planckian at low frequencies with temperature 
2h = K a /2ir until ui approaches w max , where it vanishes. 
However, in the limit |a| — > 1, K a — > and deviations 
from Planckianity appear at lower and lower frequency. 
Then, for a > 1, the flow is everywhere subsonic and 
there is no particle production because there are no neg- 
ative norm modes with uj > 0. On the contrary, when 
a < —1, the flow is everywhere supersonic. A new 
critical frequency w m ; n < w max appears, such that, for 
< uj < w m ; n , there now exist 4 asymptotic in and out 
modes, two of them with positive norm and two with 
negative norm. Particle production will thus occur with- 
out sonic horizon, but the spectrum no longer diverges 
as 1/uj for uj — > 0. 



Here we consider only the case where the horizon is 
present (|a| < 1), such that the effective temperature 
T e g is well defined for u> — > 0. The step- like analysis can 
be generalized to this situation, yielding (see Eq. (15)) 



3V3 



step 



1 



1- [q --)(!- a)D + 0(D 2 ) 



. (C4) 



where CJ max is obtained from Eq. (12) through the re- 
placement D i — ^ D(l — a). Analogously one can compute 
the effective temperature of leftgoing particles 



3^ 



step 2 

l 



,(l + a)D 2 i 



(1 - a)D + a 



2q 2 - 2q - 
2q - 1 



1 

2 

- 1 D 
~2~ 



0(D 2 ) 



(C5) 



Since the first term in the expansions of T^ ep have the 
same dependence on a as the first term of T s tep, their 
ratio 



R s t cp = 4 [q - - D 



(C6) 

docs not depend on a for small values of D. 

In Fig. 8, Teff/Tn and T e ff/T st e P are plotted for a veloc- 
ity profile with a = —0.5. Nothing qualitatively changes 
in these plots with respect to the simpler case with a = 0. 
Even if a is relevant in determining the temperature of 
the emitted spectrum both in the Hawking and the step- 
like regime, the transition between them is still governed 
by w ma i exactly as for a = 0. Thus it is again insensitive 
to the value of A/k. We have extended our study for val- 
ues of a closer than —1, and we observed no significant 
change. 
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log(A//t) log(A/0 



FIG. 8: T e s/Tn (left panel) and T c ff/r st0 p (right panel) as functions of A/o; max and oj max /K evaluated in the asymmetric flow 
profile of Eq. (CI) for a = —0.5. The validity domains of both the Hawking and the step-like regime are not affected by the 
introduction of the parameter a which governs the asymmetry of the sub and supersonic regions. 
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